Anderson impurity in pseudo-gap Fermi Systems 
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We use the numerical renormalization group method to study an Anderson impurity in a conduction 
band with the density of states varying as piui) oc with r > 0. We find two différent fixed 
points: a local-moment fixed point with the impurity effectively decoupled from the band and 
a strong-coupling fixed point with a partially screened impurity spin. The spécifie heat and the 
spin-susceptibility show powerlaw behaviour with différent exponents in strong-coupling and local- 
moment régime. We also calculate the impurity spectral function which diverges (vanishes) with 
\ulj\~'^ (I'^D in the strong-coupling (local moment) régime. 

PACS 75.20.Hr 



I. INTRODUCTION 

The behaviour of a magnetic impurities in metals is one 
of the best studied problem in condensed matter theory 
101 . In most cases it is a very good approximation to 
replace the conduction density of states by a constant as 
small variations of the density of states do not lead to 
a qualitative change of the physical properties (like the 
complète screening of the inpurity spin by the conduction 
électrons). 

The question, whether thèse physical properties are 
différent when the impurity is coupled to a Fermi System 
with a power-law density of states p{uj) ex \ujY near the 
Fermi-level was first discussed by Withoff and Fradkin 

i- 

A number of Systems are expected to show this pseu- 
dogap density of states. Among thèse are certain heavy- 
fermion superconductors where the exponent r can 
take the values r = 1 or r = 2 depcnding on the symme- 
try of the gap function. Othcr candidates are semicon- 
ductors whose valence and conduction bands touch at the 
Fermi level Q. In quasi one-dimensional metals, which 
can be viewed as realizations of the Luttinger model, the 
exponent r is a function of the Coulomb interaction |^] 
and can take values between r « 1 and r > 1 . 

Recently, the numerical renormalization group method 
(NRG) has been applied by Chen and Jayaprakash 
|g| (referred to as CY) and Ingersent to the model of an 
impurity spin coupled to a conduction band with a power- 
law density of states. In principal, this Kondo model can 
be related to a corresponding Anderson model in the limit 
of J ^ via a standard Schrieffer-Wolff transformation 
i pÔf . There is, however, a transition between a strong- 
coupling (SC) fixed point and a local-moment (LM) fixed 
point for the Kondo model at finite J so that it is a priori 
not clear whether the behaviour at this transition will be 
the same in the Anderson version of the model. 

The results of CY and Ingersent can be summarized 



as follows. For any J > Je the System approaches some 
kind of SC fixed point with the différence to the stan- 
dard Kondo model (?' — 0) that the impurity spin is not 
completely screened (a residual magnetic moment of r/8 
always remains in the zero-temperature limit). This can 
be qualitatively understood from the gradually decreas- 
ing density of states of the conduction électrons at the 
Fermi level which are responsiblc for the screening. 

The thcrmodynamic quantifies show non-Fermi liquid 
behaviour in the SC régime 



lyT) = ex T 

Xs{T) = ^T-i 



aT- 



bT- 



(1) 
(2) 



(with a,b = const.). The critical line Jc{r) starts lin- 
early for small r but diverges at r = i. In addition, 
Ingersent has shown that this divergence only holds in 
the particle-hole symmetric case and that a finite Je is 
restored away from this symmetry. (This réduction of Je 
has implications for the observability of the crossover in 
expérimental situations.) 

For any J < Je, the System approaches the LM fixed 
point where the impurity is effectively decoupled from 
the conduction band and a residual magnetic moment of 
1/4 remains. The thermodynamics in this régime have 
not yet been investigated. 

In this paper, we want to study the behaviour of an An- 
derson impurity in a pseudo-gap fermion System where 
we restrict ourselves to the symmetric case. In Sec. II, we 
want to describe our approach to the generalization of the 
NRG with a non-constant density of states, and outline 
the différences to that of CY and Ingersent. The result- 
ing formula for the hopping matrix éléments of the semi- 
infinite chain for ail n is given in Sec. III. The numerical 
results for static properties and the spectral function are 
discussed in Sec. IV and V, respectively. 
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II. GENERALIZATION OF NRG TO 
NONCONSTANT DENSITY OF STATES 

The Hamiltonian we want to study in this paper is the 
conventional single-impurity Anderson model 

H = Y.^if-iJ-io + C//1it/-iî/1u/-U 

CT 

+ E ^fc4.C/ca + E ^(Sfe) (/-laC/c. + ctZ-K^) • (3) 
ka ka 

In the model (^), é^l dénote standard annihilation (cré- 
ation) operators for band states with spin a and energy 
£fc, /l^i CT those for impurity states with spin a and energy 
£f. The Coulomb interaction for two électrons at the im- 
purity site is given by U and both subsystems are coupled 
via an energy dependend hybridization V{ek) |ÏT1 . 

In the foUowing we show that the Hamiltonian (||) is 
équivalent to a form which is more convenient for the 
dérivation of the NRG équations 

^ = E ^f/-i./-i- + C^/1iî/-it/Îu/-U 
+ E / d£g(£)4CTaeCT 

^ — 1 



(4) 



where we introduced a one-dimensional energy représen- 
tation for the conduction band with band cut-offs at ±1, 
dispersion ^(e) and hybridization h{e). The band op- 
erators fulfil the standard fermionic commutation rules 

[al^,ae'cr'] ^ S{e - e')ôaa'- 

To establish the équivalence of the Hamiltonians (||) 
and (^) we prove that for a spécifie choice of g{s) and 
h{e) they lead to the same effective action for the impu- 
rity degree of freedom. This effective action is obtained 
by integrating over the conduction électron degrees of 
freedom. For the Hamiltonian (^) one gets 

Çj EV^L+lV'am-lE^(^fc)'Gn™(fc), (5) 



(see for example |12 ). n and m count the steps on the 
imaginary time axis [0, P] with N the number of steps. ip 
and '0^ are Grassmann numbers corresponding to the im- 
purity operators. describes the unhybridized impurity. 
The G^„(fc) are Green functions for the free conduction 
électron System. 

The action corresponding to the Hamiltonian (|^) can 
be written as 



E / dexLn((l - ;^5(e))xe<Tn-i - Xe 



"§E / de/l(£) xÎCTwV'an-l + V'LXe<T«-l ■ (6) 



TV ■ 



x\an and Xean are Grassmann numbers corresponding to 
the conduction électron operators and o^ct • Integrat- 
ing over the conduction électron degrees of freedom leads 
to 



Agis)). 



(7) 



To compare the effective actions (|]) and the sum over 
k in (|^) has to be transformed to the energy intégral 

Y,V{ekfGl„Ak) = f d£F(£)V(£)G™„(£)- (8) 

This also defines the density of states for the free con- 
duction électrons p(£). The équivalence of the effective 
actions and (Q) leads to the condition 

'^d5^/i(£(5))'G^„(5) ^ fàeV{efp{e)Gl^{e). 

(9) 



This can only be fulfiUed for 

de{x) 



dx 



h{e{x)Y = V{xYp{x), 



(10) 



(with £(a;) the inverse of g{e)). For a given A(x) = 
■nV (xY p{x) there are obviously many ways of dividing 
the energy dependence between e{x) and the dispersion 
h{£{x)). One possibility is to choose 



?(£)=£ and /i(£)2 = -A(£). 



(11) 



For A(£) — A eq. ( pl] ) corresponds to the standard case 
(see eq. (2.4) in [Q). It might also be convenient to set 
h{e) ~ h. Together with the condition £(—1) = — 1 and 
£(1) = 1 this leads to 

1 [3 

eig) = — 1 H — r / da;A(a;) and 

TT/l- 



^ /' d£A(£) 

J-1 



(12) 



This équations also reduce to e{g) = g and = ^A 
for a constant A(£) = A. Equations ( pT| ) and ( p^ ) have 
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already been derived by CY ||ï^ . In a subséquent publi- 
cation p thèse authors use eq. (|Ï2|) for the mapping of 
the Kondo model on a semi-infinite chain (see Appendix 
A for a discussion of the resulting hopping matrix élé- 
ments) . 

The first possibility eq. ( pl] ) has a conceptual disad- 
vantage arising from the logarithmic discretization of the 
conduction band. Within each interval [xn+i,Xn\ and 
[—Xn, ~Xn+i], with Xn — A~", the conduction électron 
operators are expressed in terms of a Fourier expansion. 
As long as h{e)'^ is constant in each interval, the impu- 
rity couples only to the average component {p = 0) of the 
conduction électrons. Therefore it is reasonable to ne- 
glect ail the p^O-states (this becomes exact in the limit 
A — > 1). This line of reasoning obviously does not hold 
for eq. (|n|). 

On the other hand, the energy dependence of A(e) can 
be taken into account in the hybridization by defining 
h{e)'^ as the mean value 



de; 



de, 



de; 



de, 



(13) 



(14) 



(with dn = Xn — Xn+i) in each interval of the logarithmic 
discretization. This is so far not an approximation as 
the remaining energy dependence will be incorporated 
in the dispersion. The advantage of an energy dépendent 
hybridization as in eq. ( |Ï3| ) is that the resulting dispersion 
has the form g(±Xn) = for ail n, i.e. at ail points 
Xn of the logarithmic discretization. This " linear" form 
(for intermediate values g{e) = e is not fulfilled) leads to a 
scaling behaviour of the hopping matrix cléments (see eq. 
(p7|)) of the form tn oc A~"/^, slightly modified due to the 
structure of A(e). The représentation eq. (12), however, 
leads to a scaling with an effective Acfi not equal to A 
which might even dépend on the number of itérations 
thus making the analyses (of the fixed points, the relevant 
energy scale, etc.) more difHcult. 

For thèse reasons, we take the représentation eq. ( p^ ) 
in the foUowing. This gives for the hybridization part of 
the discretized Hamiltonian 



-ffhyb 



/-Icr/Oo- + /ocr/-lo- 



with 



foa = -7^ y] [in^na + 7n ^na] 



(15) 



(16) 



±\2 



(7^) 



y deA(e). 



(18) 



The discrète conduction électron operators a„(j (bna) for 
positive (négative) e correspond to those introduced in 
[^,0. According to the differential équation (|ï^) we 
should now have to solve for e{x) and invert e{x) to ob- 
tain the dispersion x{e) = g{e). This is actually not 
necessary because the single-particle énergies in the con- 
duction électron part of the discretized Hamiltonian 



na 

only dépend on the intégral over g{e) 

r± 



1 

dn 



deg(e). 



(19) 



(20) 



It can be shown that the discrète énergies are given 
by 



/^deA(e)e 
/^deA(e) ' 



(21) 



This équation, together with the form of the hybridiza- 
tion part has already been used by Sakai et al. [|^, al- 
though no dérivation was given in their article. 

The discretized Hamiltonian for the singie-impurity 
Anderson model now takes the form 

H = + t^./!iT/-iî/!u/-U 



[/llo-/o<T + îlaî-1 



III. PSEUDOGAP DENSITY OF STATES - 
MAPPING ON SEMI-INFINITE CHAIN 

We now consider a A(ù;) of the form 

A(tj) = Ao|wr, -l<w<l. 



(22) 



(23) 



The discrète énergies of the conduction électrons and 
the hybridization matrix éléments 7^ between impurity 
and the conduction électrons take the form 



21- A-C^'+i)^ 



(24) 



and 



Co = E i^^f + (7« )^) - f d.A(e), (17) (7,1)^ - {,nf ^ ^A-"(-^) (l - A-(-^)) . (25) 
», — 1 
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The mapping of the discretized Hamiltonian (^2|) onto 
the semi-infinite chain form 

^ = E ^f/-i./-i- + t^/liT/-iT/lu/-U 



OC 

an—O 



fnafn+la + fn+lafncr 



— \f-lcrfoa + foaf-lo 
TT L 



(26) 



(27) 



(Co = 1^) is described in Q and 0. The only différence 
appearing here is the r-dependence of the and 7^. 
(Note that in the non-symmetric case additional terms 
of the form Snfnafntr are generated.) For the hopping 
matrix éléments t„ we find the foUowing expressions. 

r + 1 1 - A-(''+2) 



tr. 



X i_A-(2"+''+i 
for even n and 

t„ = A-(«+'-)/2 

_ ^-(2ri+r+l 



2 1 - A-(r+i) 

1/2 



1 _ A~("+''+^) 

l _ ^-(2n+r+3) 



1/2 



(28) 



r + 1 1 - A- (■'■+2) 

r + 2 1 - A-^-'+i) 
-1/2 



1- A-("+i 

-1/2 



2^ _ ^-(2n+r+3) 



(29) 



for odd n. The équations (^8|) and (|29|) have been verified 
numerically and by analytical calculation of Iq and ti. In 
the limit n ^ 00 (pSl) and (p9|) reduce to 



in 



+ 11- A-('-+2) 



2 1 - A-('^+i) 



A" 



-ri/2 



-r/2 



n even 
n odd 



(30) 



This limit of the hopping matrix éléments has also been 
found by Ingersent [|j although the formula for ail n is 
not given in his paper. The resuit obtained by CY is 
discusscd in the appendix. 

An analytical form of the t„ for ail n > can only be 
given when the powerlaw A(ti;) = Ao|ci^|'^ extends to the 
band edges. In any expérimental realization, however, 
we expect this powerlaw only to be valid near the Fermi 
level. On the other hand, numerical studies show that 
any déviation from the form (23) close to the band edges 
merely affects the first coefficients, while the asymptotic 
behaviour again dépends on r only and is given by (30). 
Thus the qualitative behaviour near the possible low tem- 
pérature fixed points is not affected by the exact form of 
A(cj) away from the Fermi level. 



IV. RESULTS FOR STATIC PROPERTIES 



The Hamiltonian ( pq ) is solved with the NRG for the 
parameters £f = —U/2 = 10~^, A = 2.5 and différent val- 
ues for r and Aq. At each itération step we keep « 500 



States which is sufRcient for the calculation of thermody- 
namic properties. 

We first want to discuss the phase-diagram of Fig. 1 
where we have plotted the critical value A,, versus r. 
For any Aq > Ac the System flows to a strong-coupling 
fixed point (SC) similar to the fixed point in the standard 
case [Q. The energy spectrum at this fixed point can be 
explained by removing the first conduction électron site 
from the chain due to its strong coupling to the impurity. 
The remaining chain, however, has a différent structure 
as compared to the r = 0-case. Therefore this SC fixed 
point has not the Fermi liquid properties of the stan- 
dard single-impurity Anderson model (see below). For 
Aq < Ac the System always flows to the local-moment 
fixed point (LM) with the impurity effectively decoupled 
from the conduction band. Again, the resulting energy 
levels are in agreement with those of the free conduction 
électron chain. 

For both the Kondo model and the Anderson model 
Ac(7') diverges at r = ^ and we find for the Anderson 
model a logarithmic divergence 



Ac,a(?') oc 



In ( r 

2 



(31) 



However, the behaviour of Ac(r) for < r < i is quite 
différent for both models. Ingersent finds an extended 
linear région Ac(r) ex r which is approximately valid up 
to values of r = 0.4 (see inset of Fig. 1). In our case, 
Ac(r) also starts linearly and is in agreement with the 
resuit for the Kondo model up to r w 0.02, but increases 
far more rapidly for larger r. 

The différence to Q is mainly due to the fact that for 
the parameters used herc, the f-level lies within the pseu- 
dogap density of states. Under the assumption that the 
relevant coupling A' for this problem is (approximately) 
the value A(w = ef ) we have the exponential dependence 

Ao|efr = Aoe'-'"l-'^l. 



A'(r) 



(32) 



As ln|ef| has a large négative value, A'(r) is strongly 
supressed for increasing r so that a much larger Ag is 
needed to reach the strong coupling fixed point. This 
increase of the parameter régime in which local moment 
formation is observed has also been found by Gonzalez- 
Buxton and Ingersent lïsl who applied a poor man's scal- 
ing approach to the Anderson version of the pseudogap 
problem. To show that eq. (|3|) basically explains the 
différence between the Kondo model and the Anderson 
model, we have plotted in the inset of Fig. 1 both Ac,k('') 
for the Kondo model and A^^(r) — Ac,A(f )-exp(— 7.9-r). 
(the value 7.9 was chosen in order to fit A^^(r') to 
^c,k('')) The linear région of A^^(r) now extends to 
r«0.4. 

The remaining différence between Ac,k(?') and Ac.a(^) 
is due to the fact that the Kondo model and the Ander- 
son model are related via the Schrieffer-Wolff transfor- 
mation |10| only in the limit J —^ (corresponding to 
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V'^ /U 0). Therefore, the agreement of the results for 
both models is only guaranteed for A ^ 0. Away from 
the line A = 0, there is no exact mapping between the 
Kondo version and the Anderson version. 

The critical coupUng Ac(r) is determined as foUows. 
Fig. 2 shows the température dependence of the effective 
magnetic moment for U — 0.001, £f = —U/2, r = 0.48 
and différent values of A. In this graph, the LM fixed 
point (characterized by //j-cs = ^J■cttiT 0) = 1/4) is 
reached within the given température range for A — 0.01 
and A = 0.02. The value fires = r/8 = 0.06 corre- 
sponding to the SC fixed point is clearly approached for 
A = 0.16 while for A = 0.04 this value should be reached 
at a much lower température. From Fig. 2 we détermine 
Ac(ï' — 0.48) as ~ 0.03 and repeat this procédure for dif- 
férent values of r. Similar results for /icfî(T) have been 
obtained by Ingersent and CY. 

The value /iros = r/8 at the SC fixed point can also be 
derived directly from the semi-infinite chain form of the 
free conduction électron chain at this fixed point. One 
simply has to compare the effective magnetic moment for 
the System with and without the first conduction électron 
site. 

The température dependence of the spécifie heat coeffi- 
cient 7(r) = C(r) /T in the SC régime is shown in Fig. 3. 
The low température behaviour of j{T) is described by 
a power law of the form 



V. RESULTS FOR THE SPECTRAL FUNCTION 



7(T) = ciT-- + C2T- 



(33) 



Although eq. ( |3^ ) resembles an expansion in T^"^, there 
cannot be any terms like T~^^, y-4r ^^^^^ corre- 
sponding entropy would then diverge for T — > (for e.g. 
r = 0.4). 

The exponent a dcfined by "f(T) oc is shown in the 
inset of Fig. 3. In an intermediate température régime, a 
approaches the value — r consistent with the resuit of CY. 
However, for lower températures another term with the 
exponent a — —2r is dominating. This term is strongly 
suppressed in the intermediate régime due to C2 « ci. 
This crossover from the to the T~^^ behaviour is 
not due to a crossover to a new low température fixed 
point. For the spin susceptibility we confirm the resuit 
given by CY: 



X(T) = -T-i + c[T-'- 4 

In the LM régime we find 

7(T) = C3T'^-i. 

and 



(34) 



(35) 



(36) 



The impurity spectral function 



X e 



J -la 



-^^") , (37) 



(with the partition function Z = J^m^^'Pi"!^-^'"'-))^ 
not yet been calculated in the previous papers on the 
pseudogap problem. We assume that the groundstate 
energy Eg is set to zéro and concentrate on the zero- 
temperature limit, in which the spectral function takes 
the form 



E 



("9 


J-la 




(rie 


J-la 


rug^ 



Ô{UJ + ErnJ 



Here, the partition function Z equals the total degener- 
acy of the groundstate. The n 



label ail states with 
energy E = Eg~0 and the rte, nie label the excited states. 



9' ' 



The first term in eq. ( |38[ ) would correspond to a transi- 
tion between différent states with E — 0, but such a term 
(resulting in a i5-function at the Fermi level) is not présent 
in the NRG results. There is one state with excitation 
energy E^x — > for iV — > oo but its matrix clément with 
the ground state vanishes as iV ^ oo. 

In order to obtain the full frequency dependence of 
the spectral function within the NRG, it is necessary to 
combine the information of ail itération steps as in each 
itération the results are only given for a certain frequency 
range (see also |^,0). 

In Fig. 4 we show results for the spectral function for 
r = 0.25, A > Ac (solidline, SC régime), r = 0.25, A < Ac 
(dotted line, LM régime) and r = 0.75 (dashed line, LM 
régime). For thèse calculations we used A = 2 and kept 
w 800 states at each itération. 

We find that A{uj) diverges as |a;|~'" for w — > for any 
set of parameters which lies in the SC régime. Note that 
this resuit suggests the conventional behaviour A{uj) ^ 
l/(7rA(ijj)) as a; — > for the SC case. Together with 
the resuit 7(T) ^ T^^ (neglecting the second term in 
(33) for the moment) one could be tempted to interprct 
thèse results within a standard Fermi liquid approach. 
Let us however emphasize that in spite of thèse results 
the System is not a Fermi liquid for any r > 0. 

This observation becomes more évident in the LM 
régime, where the behaviour of the spectral function is 
qualitatively différent. Namely, in contrast to the SC case 
we find that the spectral function vanishes as A(oj) oc \ùj\'^ 
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here. In addition, no qualitative différence, apart from 
the exponent, can be observed between the cases r> 0.5 
and r < 0.5. 



VI. SUMMARY 

To summarize, we have studied the problem of an 
Anderson impurity in a pseudo-gap Fermi System at 
particle-hole symmetry using a generalization of tfie nu- 
merical renormalization group method [^j^. We find a 
befiaviour similar to that of the corresponding Kondo 
model investigated by CY |^ and Ingersent 0. How- 
ever, the critical Une Ac separating the strong-coupling 
and local-moment régimes of thèse two models shows a 
quite différent form between r = (where Ac starts lin- 
early) and r — 1/2 (where Ac diverges). This différence 
is mainly due to the fact that we have chosen the f-level 
to lie within the pseudogap. In both strong-coupling 
and local-moment régime the thermodynamic quantifies 
spécifie heat and spin-susceptibility show powerlaw be- 
haviour. 

We also presented the first calculations of the impurity 
spectral function for this model. We find A{uo) ex in 
the strong-coupling régime and A{lo) oc \uj\^ in the local- 
moment régime. We do not find any indication that the 
local moment fixed points for r < 1/2 and r > 1/2 arc 
différent. 

As shown by Ingersent for the Kondo model, the crit- 
ical values Ac take finite values as soon as particle-hole 
symmetry is violated. It is of course interesting to see, 
whether this réduction of the critical coupling is the 
same for the Anderson model (work on this problem is 
in progress). 

Another interesting question is the relevance of the 
model studied here in the context of the dynamical mean 
field theory (for récent reviews see |l8yi9| ]). The effective 
single impurity Anderson model appearing in the dynam- 
ical mean field theory is coupled to a (self-consistently 
determined) effective médium. There is a possibility 
that the density of states corresponding to this effec- 
tive médium develops a pseudo-gap structure under cer- 
tain conditions (e.g. near the métal- insulator transition). 
Also, the density of states of the infinité dimensional gen- 
eralization of the honeycomb (c?= 2) and diamond (d — 3) 
lattices is proportional to IcijI near the Fermi level. 

We wish to thank J. Keller and G. M. Zhang for a 
numbcr of stimulating discussions. Onc of us (R.B.) was 
supported by a grant from the Deutsche Forschungsge- 
meinschaft, grant No. Bu965-1/1. 



APPENDIX A: OTHER DISCRETIZATION 
SCHEMES 

In this appendix, we want to show that the discretiza- 
tion used by CY leads to the same hopping matrix élé- 
ments tn apart from a redéfinition of the discretization 
parameter A. This équivalence, however, is restricted to 
the spécial form of A(w) and is not valid in gênerai. 

For A(w) = Aol^r eq. (|ï|) leads to 



aie) 



£'■+1, 
Aq 1 

n r + 1 



(Al) 
(A2) 



The discrète énergies take the form 

1 



dn 



deg{e) 



r + 1 1 - A '■+! 
r + 2 1-A-i 
r + 11- Â'i-+^) 



■2 1- A-('-+i) 



(A3) 



In the last équation we have defi ned A = A-'+i . 

The hybridization in eq. (A2) is independent of fre- 
quency, therefore the resuit for (7,^) is the same as for 
a constant A{oj) — Ao/{r + 1) 



Ao 



-A-" (1- A-i) 



With A = A('~+i) eq. (0) gives 

„±\2 _ Ao ^^„(r+l) 



(A4) 



(A5) 



Eqs. (|A3|) and (|A5D correspond to Eqs. (gj) and (|2^) 
with A bcing replaced by A. Obviously, also the resulting 
matrix cléments i„ will have exactly the same form as in 
eqs. (|8|) and ( p9| ) and in terms of A we find in the limit 
of large N 



r 2 1 - A-1 1 A 



For the ratio of i„/t„_i we then find 

tn n^ao J A^t'^+i) : n eveu 
t~^ ' I A-1/2 : nodd ' 



n even 
n odd 



(A6) 



(A7) 



corresponding to the resuit given in CY. 
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FIG. 1. r-dependence of the critical coupling Ac which 
seperates the strong-coupling régime (A > Ac) from the lo- 
cal-moment régime (A > Ac). The fiUed squares show the 
resuit for Ac,A(r). In the inset we compare the scaled criti- 
cal coupling Ac,a('") (open squares) with the resuit obtained 
by Ingersent for the Kondo version of the Hamiltonian (3) 
(circles). 
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FIG. 2. Température dependence of the effective magnetic 
moment for U — 0.001, Sf — ~U/2, r — 0.48 and varions val- 
ues of A. For A < 0.03, the System flows to the local-moment 
fixed point with the corresponding effective magnetic moment 
fj,o — 1/4. For A > 0.03 the magnetic moment is partially 
screened and approaches the residual value ^rcs — r/8 — 0.06 
from above. 
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FIG. 3. Température dcpenderice of the spécifie heat coef- 
ficient 7(r) = C{T)/T for U = 0.001, et = -U/2, r = 0.48 
and various values of A. The inset shows the température 
dependence of the exponent a{T) defined by 'y{T) oc . In 
an intermediato régime, a approaches —r as expected from 
eq. (1) but for lower températures a stronger divergent term 
with a « —0.93 « — 2r dominâtes. This behaviour is not due 
to the crossover to a new fixed point. 
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FIG. 4. Impurity spectral function for U = 0.001, 
£f = —U/2, and différent values of r and A. for r — 0.25 
and A = 0.02 (solid line, SC-regime), the spectral function 
diverges as A{u!) oc \u!\~^ for |w| — »■ 0. In the LM-regime 
(for both r = 0.25 and A = 0.0002 (dotted line) and r = 0.75 
and A — 0.02 (dashed line)) the spectral function vanishes 
as The inset shows the coefficient a{u!) defined by 

A{lo) oc 
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